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1. A quantitative analysis of multiparticle data requires using generators 
of computer "events". The trouble is that the models employed in this con- 
text are essentially classical and that their predictions are inadequate and 
sometimes turn out to be misleading in the study of those aspects of multi- 
particle production where quantum interference is important [||]. This is a 
serious problem because, as evidenced by the studies of the so-called inter- 
mittency phenomenon , the momentum-space short-range correlations are 
to a large extent dominated by the Hanbury-Brown-Twiss (HBT) effect [0], 
i.e. by quantum interference. 

We would like to emphasize that this does not mean at all that these 
correlations are "trivial" as is sometimes claimed. On the contrary, since 
quantum interference is sensitive to the space-time development of the colli- 
sion process it yields precious information thereon and its effects should be 
examined with most attention. 

The aim of this communication is to outline a way of implementing the 
effects of quantum statistics in Monte-Carlo simulations of multiparticle phe- 
nomena. One could wonder whether this is possible at all, since the Monte- 
Carlo method deals with probabilities while quantum interference is a conse- 
quence of adding amplitudes whose phases are essential. We will show, how- 
ever, that one can devise a systematic and practical approach to the problem 
if Wigner functions and not scattering amplitudes are used to describe the 
multiparticle system. 

2. We shall now write the n-particle spectrum in terms of the Wigner 
function, first for distinguishable and then for identical secondaries. 

Let ip{q,; a) denote a wave- function describing a stationary state produced 
at high-energy. Here q refers to the momenta qi,q2, ■■.,qn of the n produced 
spinless particles and a denotes all other parameters, assumed irrelevant for 
our problem. For the moment, the secondaries are supposed to be distin- 
guishable. 

The n-particle spectrum is, of course 

^o(q) = I] I ^(q;a) q = (gi, (i) 

a 

The weights generated by a Monte-Carlo algorithm are directly proportional 
to fio(q). 

Going over to the coordinate representation we write 
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f]o(q) = / (ix(ixVq-(^-^') p(x; x'), x = {x,, Xn) (2) 

where 

x') = J2 a);A*(x'; a) (3) 

a 

is a density matrix, ip being the Fourier transform of i/j. Eq. can further 
be rewritten as 

Qo{q) = y'dx+Vl/(q;x+) (4) 
where [notation: x"*" = |(x + x'), x^ = x - x'] 

iy(q;x+) = J (ix-e*^-^Xx; x') (5) 

is the generahzed Wigner function [Q, the quantum analog of the classical 
Boltzmann phase-space density. It is real and it gives the observable spec- 
trum when integrated, as seen in (^. 

Let us now assume that the secondaries are all identical (the generalisa- 
tion of the discussion to the case where there are several species of identical 
secondaries is straightforward). Let V denote some arbitrary permutations 
of the integers 1, n and q^ the corresponding permutation of the momenta 
gi, Once the wave function ?/^(q; a) is symmetrized with respect to the 
momenta of produced bosons, eq. (|) becomes 

Q{q) = ^ E / rfx(ix'e'(^-q^-^'-q^') p(x; x') (6) 

which is rewritten as 

^(q) = ^ E / rfx+e^^" (q--q-') T4^(51^±^; x+) (7) 

with VF(q;x"'") defined in (^). Thus we observe that the same Wigner func- 
tion determines the spectrum before and after the symmetrization. This has 
a simple physical reason: All information, compatible with the rules of quan- 
tum mechanics, about what happens in the full phase-space of coordinates 
and momenta (in contrast to the momentum-space alone) is encoded in the 
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Wigner function. This information includes the phases of different waves and 
therefore all that is needed to predict the interference patterns. A further 
advantage of the Wigner function is that it appeals directly to one's intuition. 
A few words of caution are necessary at this point, however. 

Contrary to the Boltzmann phase-space density, the Wigner function is 
locally not positive definite. It can actually oscillate quite violently. The 
oscillations integrate to zero in (|]) but can conspire with oscillating terms 
in the integrand of (|^) to contribute significantly to the result. This is how 
quantum mechanics shows up in the problem. Thus, it is clear that regard- 
ing the Wigner function as a phase-space density is possible only when the 
function is smoothed by averaging the oscillations out. Physically it means 
appropriate smearing of coordinates and momenta. The price to pay is that, 
in general, the resulting probabilistic model can only be trusted when the 
momentum differences appearing in (|^) are not too large. 

3. The standard Monte-Carlo algorithms rest on models of momentum- 
space densities. The goal to achieve is to correct the weights of Monte-Carlo 
events, once they have been generated according to the distribution fiolq)- 
Our proposal consists in going from f2o(c[) to f2(q) by modelling the Wigner 
function. 



we see that, if the Wigner function is regarded as a phase-space density, 
w(q; x+) has the meaning of a conditional probability: given that the parti- 
cles with momenta qi, are present in the final state, w is the probability 
that they are produced at the points xf , ...,x^. The problem is to construct 
a viable model for w{q;, x^). 

In the absence of additional information it seems reasonable to start with 
the working assumption that the likelihood to radiate a particle from a given 
space point is statistically independent of what happens to other particles. 
This means that w factorizes: 



Writing 




(8) 




(9) 



where 




(10) 
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Substituting (D-(D into (g) one finds 



^(q) = ^ E ^o(^^^) n^fe' i<lv - <lv')j] (11) 
n. J, J,, / J 

where 

w{q, A) = / d^x e^-'^wiq, x) (12) 

Eq. (p!OD and the reahty of w{q, x) imply w{q, 0) = 1 and w{q, A) = w*{q, —A), 
respectively. This guarantees that fi(q) calculated from (|TI|) is real. 

Eq. ( p!TD can be used as it stands when one has an explicit formula for 
Oo(ci)- In practice, however, fio(ct) is constructed iteratively by a Monte- 
Carlo algorithm and at a given stage of the simulation it is computed for 
one configuration of momenta. To deal with this complication we observe 
that one does not make a big error by replacing in (|TI]) ^^o( ^^"^^^' ) by 
f2o(q-p). Indeed, those terms in eq. ([TI|) where this approximation is poor 
are suppressed by the rapidly decreasing factors w and thus need not be 
calculated with a great precision. The equation (|1T]) now becomes 

^(q) = 3E{^o(qp) E n^fe' (Qv - qv')j]} (13) 

-p -pr j 

Thus, once a configuration of momenta has been generated by the original 
algorithm, the weight of the event in question has to be multiplied by the 
correction factor Re{J2v' Yij ^['Zji ilv — Qp')j]} in order to take care of the 
HBT interference. This is the result sought. In ( |TB| ) the sum over V just 
expresses formally the fact that even in a classical model the labelling of 
identical particles is merely a matter of convention. 

The function w{q, A) is unknown. It can either be taken from a model 
or, perhaps more reliably, be determined by fitting 2-body HBT correlations. 
The weight of an event is then completely determined and the procedure can 
be used to study the implications of quantum statistics for other aspects of 
the production process. 

4. Adopting the probabilistic interpretation of the Wigner function, the 
proposed approach allows an intuitive interpretation of the results in terms 
of the space-time structure of the region of particle emission. As already 
mentioned this is meaningful when momentum differences are not too large. 
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It should be clear that a simple physical meaning can be ascribed to 
w{q, x) rather than to w{q, A). Let us briefly outline how one might proceed 
in modelling w{q, x). We imagine that a particle with momentum q is emitted 
from a diffuse source centered at X — Xq (q). We are thus led to write 

w{q,x) = P[q,x - Xo{q)] (14) 
The simplest choice for P{q, x) would be to take a Gaussian 

Pg{x) = (15) 

TTiaxayaz 

with a = criq), to take into account the possible dependence of the shape of 
the source on q. More generally one can set 

Piq,x) = J d^a H{q,a)PGix) (16) 

This gives 

w{q, A) = e^^-'f^") J d^'a H{q, a) e-i(-'^'+-?^?+-'^'), (17) 

a formula that seems general enough to accomodate all physically reasonable 
choices for w. 

Notice that, for obvious physical reasons, the distribution ([T^) can only 
depend on differences xo{q) — xo^q'). This is indeed the case as can be seen 
by observing that Y.jilv ~ 1v')j = for any two permutations V and V. 

A further simplification is obtained by assuming that the place where a 
particle is produced depends at most on some global characteristics of the 
collision, like the total energy, to give an example. This is presumably a good 
assumption as long as the source can be regarded as static. Technically, it 
amounts to neglect the first argument of w{q, A). 

5. A few comments are in order: 

(i) The proposed approach, even used in its simplest version, enables one 
to incorporate into a Monte-Carlo simulation the collective nature of the HBT 
effect. 

(ii) Our ansatz is to be considered as a working assumption to be 
upgraded when more information is available (e.g. if two identical particles 
result from the decay of the same resonance). 
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(iii) The ansatz can be checked against the data on higher order correla- 
tions. This may lead to a discovery of hitherto entirely unknown correlation 
structures in the space-time development of the collision and is, therefore, 
potentially of great interest. 

(iv) The motivation of our probabilistic approach to the Wigner function 
is mostly phcnomcnological. It would be, of course, highly desirable to invest 
more effort in developing theoretically based models of the Wigner function. 

To summarize, we have shown that the effects of quantum interference in 
multiparticle production can be naturally expressed in terms of generalized 
Wigner functions. Such formulation allows one to incorporate these effets 
into Monte-Carlo generators and to give them an intuitive interpretation in 
terms of the space-time development of the interaction. It also explicitly 
demonstrates that in order to obtain theoretically founded predictions for 
the HBT correlations one should compute the Wigner functions from the 
underlying theory. 
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